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Abstract: We start by presenting the current status of a complex flavour conserving 
two-Higgs doublet model. We will focus on some very interesting scenarios where unex¬ 
pectedly the light Higgs couplings to leptons and to b-quarks can have a large pseudoscalar 
component with a vanishing scalar component. Predictions for the allowed parameter space 
at end of the next run with a total collected luminosity of 300/6 _1 and 3000/6 _1 are also 
discussed. These scenarios are not excluded by present data and most probably will survive 
the next LHC run. However, a measurement of the mixing angle </> T , between the scalar 
and pseudoscalar component of the 125 GeV Higgs, in the decay h -A r + r~ will be able to 
probe many of these scenarios, even with low luminosity. Similarly, a measurement of <4>t in 
the vertex tth could help to constrain the low tan /? region in the Type I model. 
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1 Introduction 

The discovery of the Higgs boson at the Large Hadron Collider (LHC) by the ATLAS [1] 
and CMS [2] collaborations has ignited a very large number of studies in the context of 
multi-Higgs models. It is now clear that some features of the Higgs couplings to fermions 
and gauge bosons have to be well within the Standard Model (SM) predictions. Also, even 
if other heavy scalars are far from being experimentally excluded, there is still no hint of 
scalar particles other than the 125 GeV one. However, even if no large deviations from 
the SM were found, many of its extension are still in agreement with all experiment data. 
Many models provide interesting scenarios that can be probed at the next LHC run while 
contributing to solve some of the outstanding problems in particle physics. Such is the 
case of the complex two-Higgs double model (C2HDM). The 2HDM was first proposed by 
T. D. Lee [3] as an attempt to understand the matter-antimatter asymmetry of the universe 
(the 2HDM is described in detail in [4, 5]). 

The 2HDM is a simple extension of the SM where the potential is still invariant under 
SU(2) x U( 1) but is now built with two complex scalar doublets. The complex two-Higgs 
doublet model is the version of the model that allows for CP-violation in the potential, 
providing therefore an extra source of CP-violation to the theory. The existing experimental 
data and in particular the one recently analysed at the LHC has been used in several studies 
with the goal of constraining the parameter space of the C2HDM [6-9] or just the Yukawa 
couplings [10]. 

The main purpose of this work is to analyse scenarios in the C2HDM that deviate from 
the SM predictions, while being in agreement with all available experimental and theoretical 
constraints. These are scenarios where the scalar component of the Higgs coupling to 
leptons or to b-quarks vanishes. The respective pseudoscalar component has to be non-zero 
which does not necessarily imply a very large CP-violating parameter. Even if the scalar 
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component is not exactly zero, there are still Yukawa couplings where the pseudoscalar 
component can be much larger than the corresponding scalar component. 

We will start by discussing the status of the C2HDM. Presently the processes pp —»> 
h WW(ZZ ), pp ^ h ^ 77 and pp —» h —» r + r _ are measured with an accuracy of 
about 20%. On the other hand pp —)• V(/i —)► 66 ) has been measured at the Tevatron and at 
the LHC with an accuracy of about 50% [11, 12 ] while for pp h Z 7 an upper bound 
of the order of ten times the SM expectation at the 95% confidence level was found [13, 14]. 
In order to understand how the model will perform at the end of the next LHC run we 
use the expected precisions on the signal strengths of different Higgs decay modes by the 
ATLAS [15] and CMS [16] collaborations (see also [17]) for = 14 TeV and for 300 and 
3000 /6 _1 of integrated luminosities. As previously shown in [9], the final states VV, 77 
and t + t~ are enough to reproduce quantitatively the effect of all possible final states in the 
Higgs decay. Therefore, taking into account the predicted precision for the signal strength, 
we will consider the situations where, at 13 TeV, the rates are measured within either 10% 
or 5% of the SM prediction. We should note that no difference can be seen in the plots 
when the energy is changed from 13 to 14 TeV as discussed in [9]. 

This paper is organized as follows. In Section 2, we describe the complex 2 HDM and 
the constraints imposed by theoretical and phenomenological considerations including the 
most recent LHC data. In Section 3 we discuss the present status of the model and in 
Section 4 we discuss the scenarios where the pure scalar component of the Yukawa coupling 
is allowed to vanish. Our conclusions are presented in Section 5. 

2 The complex 2HDM 

The complex 2HDM was recently reviewed in great detail in [9] (see also [ 6 , 18-23]). 
Therefore, in this section we will just briefly describe the main features of the complex two 
two-Higgs doublet with a softly broken Z 2 symmetry </>i ^ —</>2 whose scalar 

potential we write as [5] 

Vh = raiil<^| 2 + "*22l<H 2 — 77 T.f 2 01^*2 - (rnl 2 )* (j)l(j)l 

+ I ^2 1 4 + A 3 | 0 l| 2 |^ 2 | 2 + A4 (< 4 </> l) 

+■^■(d^) 2 + 7(^2^ i) 2 - ( 2 -i) 

All couplings except m\ 2 and A 5 are real due to the hermiticity of the potential. The 
complex 2HDM model as first defined in [18], is obtained by forcing arg(As) % 2arg(m^ 2 ) 
in which case the two phases cannot be removed simultaneously. From now on we will refer 
to this model as C2HDM. 

We choose a basis where the vacuum expectation values (vevs) are real. Whenever we 
refer to the CP-conserving 2 HDM, not only the vevs, but also m\ 2 and A 5 are taken real. 
Therefore, 2HDM refers to a softly broken Z 2 symmetric model where all parameters of the 
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potential and the vevs are real. Writing the scalar doublets as 


' 6l (73(^1 + ? ?i + ®xi)) ’ (73(^2 + % + %)) ’ ^ 2 ' 2 ^ 

with ?; = \Jv\ + ^2 = (v^G ^) -1 / 2 = 246 GeV, they can be transformed into the Higgs 
basis by [24, 25] 

where tan/3 = ^ 2 /^ 1 , cp = cos/3, and = sin/3. In the Higgs basis the second doublet 
does not get a vev and the Goldstone bosons are in the first doublet. 

Defining 773 as the neutral imaginary component of the H 2 doublet, the mass eigenstates 
are obtained from the three neutral states via the rotation matrix R 


/ hA 



h 2 

= R 

m 

W) 


\m) 


which will diagonalize the mass matrix of the neutral states via 

RM 2 R t = diag (m 2 , m 2 , m 2 ) , 


(2.4) 


(2.5) 


and mi < m 2 < m 3 are the masses of the neutral Higgs particles. We parametrize the 
mixing matrix R as [20] 


R = 


( cic 2 

-(C 1 S 2 S 3 + sic 3 ) 
-cis 2 c 3 + sis 3 


S\C 2 S2 ' 

C1C3 — S1S2S3 c 2«3 
-(cis 3 + S1S2C3) C 2 C 3 J 


( 2 . 6 ) 


with Si = sinc^ and q = cosc^ (i — 1, 2, 3) and 


— 7t/2 < OL\ < 7t/2, —7t/2 < 0^2 < 7T /2, — 7t/2 < «3 < 7T /2. (2.7) 


The potential of the C2HDM has 9 independent parameters and we choose as input 
parameters v : tan/3, m H ± : a±, <a 2 , ol 3 , mi, m 2 , and Re(m 2 2 ). The mass of the heavier 
neutral scalar is then given by 

2 _ mf-Ri 3 (-Bi 2 tan / g--Bn) + m% #23(^22 tan/? - i? 2 i) 

3 -^33(^31 — i? 32 tan/ 3 ) 

The parameter space will be constrained by the condition m 3 > m 2 . 

In order to perform a study on the light Higgs bosons we need the Higgs coupling to 
gauge bosons that can be written as [ 6 ] 

C = cpRn + spR \2 = cos ( 0 ^ 2 ) cos (ol\ — /3), (2.9) 

and the Higgs couplings to a pair of charged Higgs bosons [ 6 ] 

— A = C {3 [s^Ai45 + C3A3] R\i + S {3 A245 + S^X 3] R \2 + SpCp Im(As) i?!3, ( 2 . 10 ) 
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where A 145 = Ai — A 4 — Re(A 5 ) and A 245 = A 2 — A 4 — Re(A 5 ). Finally we also need the 
Yukawa couplings. In order to avoid flavour changing neutral currents (FCNC) we extend 
the symmetry to the Yukawa Lagrangian [26]. The up-type quarks couple to 02 and the 
usual four models are obtained by coupling down-type quarks and charged leptons to 02 
(Type I) or to 0i (Type II); or by coupling the down-type quarks to 0 1 and the charged 
leptons to 02 (Flipped) or finally by coupling the down-type quarks to 02 and the charged 
leptons to 0i (Lepton Specific). The Yukawa couplings can then be written, relative to the 
SM ones, as a + 2675 with the coefficients presented in table 1. 



Type I 

Type II 

Lepton 

Specific 

Flipped 

Up 

lf-^lfl5 

if-^lfls 

if-fc/llf-B 

ElZ — j r „El3^ 
s f 3 lC P s p T5 

Down 

if + **lfl6 


if + “Ulf'75 

if-^ifis 

Leptons 

if + fc 0lfl5 

if -“S if 75 

if - “8 if 75 

if +<<»lfl5 


Table 1. Yukawa couplings of the lightest scalar, /q, in the form a + 2675 . 


From the form of the rotation matrix R (2.6), it is clear that when S 2 = 0, the pseu¬ 
doscalar 773 does not contribute to the mass eigenstate h\. It is also obvious that when 
<§2 — 0 the pseudoscalar components of all Yukawa couplings vanish. Therefore, we can 
state 


1 52 1 =0 =^> h\ is a pure scalar, (2.11) 

|<§ 2 1 = 1 h\ is a pure pseudoscalar. (2.12) 

There are however other interesting scenarios that could be in principle allowed. We 
could ask ourselves if a situation where the scalar couplings ap « 0 (F = U,D,L) is still 

allowed after the 8 TeV run. As ajj is fixed (the same for all Yukawa types) and given 

by ajj = R 12 /sp = siC 2 /sp, it can only be small if s\ ~ 0. If instead C 2 ~ 0 the h{VV 
coupling C in eq. (2.9) would vanish which is already disallowed by experiment. There is 
one other coupling that could also vanish, which is Rn/cp = c^/c# (this is for example 
the expression for oq in Type II). Again this scalar part could vanish if <7 ~ 0. In either 
case, si ^ 0 or ci ^ 0 , the important point to note is that the pseudoscalar component of 
the 125 GeV Higgs is not constrained by the choice of oq because it depends only on <§ 2 - 
We will discuss these scenarios in detail in section 4. 

3 Present status of the C2HDM 

We start by briefly reviewing the status of the C2HDM after the 8 TeV run. We will 
generate points in parameter space with the following conditions: the lightest neutral scalar 
is mi = 125 GeV 1 , the angles 07^,3 all vary in the interval [—vr/2, 7 t/ 2], 1 < tan/3 < 30, 

1 The latest results on the measurement of the Higgs mass are 125.36 ± 0.37 GeV from ATLAS [27] and 
125.02 + 0.26 - 0.27 (stat) +0.14 - 0.15 (syst) GeV from CMS [28]. 
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mi < m 2 < 900 GeV and — (900 GeV ) 2 < Re(m\ 2 ) < (900 GeV) 2 . Finally, we consider 
different ranges for the charged Higgs mass because the constraints from B-physics, and in 
particular the ones from b —)> sy, affect differently Type II/F and Type I/LS. In Type II 
and F the range for the charged Higgs mass is 340 GeV < m H ± < 900 GeV due to b —» sy 
which forces m H ± > 340 GeV almost independently of tan/3 [29]. In Type I and LS the 
range is 100 GeV < m H ± < 900 GeV because the constraint from b —sy is not as strong. 
The remaining constraints from B-physics [30, 31] and from the R 5 = Y{Z bb)/T(Z 
hadrons) [32] measurement have a similar effect on all models forcing tan /3 > 1. The choice 
of the lower bound of 100 GeV is due to LEP searches on e + e - —)► H + H~ [33] and the latest 
LHC results on pp —>> it(—> H + b) [34, 35]. Very light neutral scalars are also constrained 
by LEP results [36]. 

All points comply to the following theoretical constraints: the potential is bounded 
from below [37], perturbative unitarity is enforced [38-40] and all allowed points conform 
to the oblique radiative parameters [41-43]. 

The signal strength is defined as 

hi _ cBR(hj —>» /) 

“ <7 sm BR SM (/ij —> /) 

where a is the Higgs boson production cross section and BR (hi f ) is the branching 

ratio of the hi decay into the final state /; a SM and BR SM (/i f ) are the respective 

quantities calculated in the SM. The gluon fusion cross section is calculated at NNLO using 
HIGLU [44] together with the corresponding expressions for the CP-violating model in [9]. 
SusHi [45] at NNLO is used for calculating bb —>> /i, while Vh (associated production), tth 
and VV h (vector boson fusion) can be found in [46]. As previously discussed we will 
consider the rates for the processes /iyy, /i 77 and p rr to be within 20 % of the expected 
SM value, which at present roughly matches the average precision at la. It was shown 
in [9] that taking into account other processes with the present attained precision has no 
significant impact in the results. 


(3.1) 


60 

30 

0 

-30 

-60 



«i (deg) 


«i (deg) 


Figure 1 . 02 vs. ol\ for Type I (left) and Type II (right). The rates are taken to be within 20% 
of the SM predictions. The colours are superimposed with cyan/light-grey for pvv, blue/black for 
(i TT and finally red/dark-grey for with a center of mass energy of 8 TeV. 
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We start by examining the parameter space for a center of mass energy of y/s — 8 
TeV corresponding to the end of the first LHC run. The rates are taken at 20% and the 
effect of considering each of the rates at a time is shown by superimposing the colours, 
cyan/light-grey (/iyy), blue/black (/i rr ) and finally red/dark-grey (/i 77 ). In figure 1 we 
present the allowed space for the angles a 2 vs. oq for Type I (left) and Type II (right) 
with all theoretical and collider constraints taken into account. The corresponding plots 
for the Flipped (Lepton specific) are very similar to the one for Type II (I) and are not 
shown. It was expected that a 2 would be centred around zero where the pseudoscalar 
component vanishes. Also a i plays the role of a + 7r/2, where a is the rotation angle in 
the CP-conserving case 2 . In previous works for the CP-conserving model [47, 48] we have 
made estimates for the allowed parameters based on the assumption that the production 
is dominated by gluon fusion and that T(h± —> bb) is to a good approximation the Higgs 
total width. Under a similar approximation, we can write for Type I and large tan/3 (when 
tan/3 > 1, ^ < 1 and we recover the CP-conserving Yukawa couplings) 

Hyv ~ cos 2 a 2 cos 2 (/3 — aq). (3.2) 


Since we are considering a 20% accuracy, it is clear that neither cos<a 2 nor cos(/3 — aq) can 
be close to zero. In fact, a measurement of fiyy with a 20% (5%) accuracy and centred at 
the SM expected value implies cos 2 a 2 > 0.8(0.95) and consequently | sin 0 ^ 2 1 0.45(0.22) 

and \ot 2 1 < 27°(13°). Although the approximations captures the features, the plot does not 
reproduce the exact value of the limit, which for a 20% accuracy is slightly below 50°. 

In figure 2 we show tan/3 as a function of Rn (left), R \ 2 (middle) and i?i 3 (right). The 
upper plots are for Type I and the lower plots for Type II. Again the differences of Type 
II (I) relative to F (LS) are small and we do not show the corresponding plots. We start 
with i?i 3 which is just sina^, thus measuring the amount of CP-violation for the 125 GeV 
Higgs, that is, the magnitude of its pseudoscalar component. The allowed points are centred 
around zero where we recover a SM-like Higgs Yukawa coupling for the lightest scalar state. 
The differences between the models only occur for large tan /3, reflecting the different angle 
dependence of the couplings in the various models. We now discuss R \ 2 = sin aq cos a 2 . 
Using the same approximation for fiyy as in eq.(3.2) we can write for large tan/3 


^vv 


TD 2 

n 12 > 


(3.3) 


which means that if we take |i? 2 2 | > 0.8 then R\ 2 > 0.89 or R\ 2 < —0.89. These are 
exactly the bounds we see in the plots for Type I. Therefore, as already happened for the 
CP-conserving case it is mainly fiyy that constrains \R\ 2 \ to be close to 1 especially for 
large tan/3. Finally Rn — cosoq cos a 2 is only indirectly constrained by the bounds on 
aq and a 2 . Since the pure scalar part of the coupling relative to the SM is proportional 
to R\ x (1 + tan 2 /3) it is natural that when Rn increases, tan/3 decreases. However, the 
most important point to note is that Rn — 0 is allowed. Although Rn is never part of 
the Yukawa couplings in Type I, it appears in pure scalar couplings for down-type quarks 

2 We can choose a parametrization where the angle a± is exactly a in that limit. See for example the 
definition of the rotation matrix in [7] as compared to our equation 2.6. 
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Type I 


Type I 



Figure 2. Top: tan/3 as a function of Rn (left), R 12 (middle) and Ris for Type I. Bottom: same 
but for Type II. The rates are taken to be within 20% of the SM predictions. The colours are 
superimposed with cyan/light-grey for /iyy, blue/black for fi TT and finally red/dark-grey for /z 77 
with a center of mass energy of 8 TeV. 


or/and charged leptons in the remaining types. This in turn implies that scenarios where 
an — 0 and/or — 0 are not excluded. Models Type II, F and LS can therefore have 
a pure pseudoscalar component for some of its Yukawa couplings. This scenario will be 
discussed in detail in the next section. 



Figure 3. tan/3 as a function of sin(ai — 7 r/ 2 ) with all rates at 20% for Type I (left), Type II 
(middle) and LS (right). All angles are free to vary in their allowed range (cyan/light-grey) and we 
impose the constraint 82 <0.1 (blue/black) and 82 < 0.05 (red/dark-grey). 


In figure 3 we present tan/3 as a function of sin(oa — 7 t/2 ) with all rates at 20% for 
Type I (left), Type II (middle) and LS (right). All angles are free to vary in their allowed 
range and we present scenarios for which 82 < 0.1 and 82 < 0.05. We plot ol\ — 7 r /2 instead 
of a\ to match the usual definition for the CP-conserving model. Since we recover the 
CP-conserving h\ couplings when 82 = 0, the red/dark-grey outer layer for Type II and 
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LS has to match the bounds for the angle a in the CP conserving case which is indeed the 
case [5]. If we identify oq with <a + 7r/2, where a is the rotation angle for the CP-conserving 
scenario, we can write the coupling to gauge bosons as 

9hvv = cos (<4 9hvv ■ (3-4) 


Hence, for Type I /ayv will either give the same bound as in the CP-conserving case or 
worse as cos (c^) decreases. However, for Type II, the same approximation that lead to eq. 
(3.2) for Type I results for Type II in 


Vvv 


COS 2 OL2 cos 2 (/3 — Oil) 
tan 2 /3 


sin 2 oq cos 2 OL 2 + sin 2 OL 2 cos 2 /3 
cos 2 oq cos 2 OL 2 + sin 2 ol 2 sin 2 /3 


(3.5) 


Again, if 52 = 0 we recover the CP-conserving expression. However, it can be shown that 
larger values of 82 together with smaller values of tan/3 still fulfil the constraints on the 
rates. We conclude that in Type I, the allowed parameter space is the same as in the 
CP-conserving case while, for the remaining types and for a given aq, the upper bound on 
tan/3 is the same as in the CP-conserving case. But, now, there is no lower bound on tan/3. 


4 The zero scalar components scenarios and the LHC run 2 

In the previous section we have shown that R\\ — 0 is still allowed, which implies that 
the pure scalar components of the Yukawa couplings can be zero in some scenarios. This 
possibility arises in Type II, F and LS. In particular for Type II we can have either ap — 0 
or cll = 0 while in F (LS) only ap = 0 (a^ = 0) is possible. For definiteness let us now 
analyse the case where — 0 in Type II. Since ajj — R\i/cp = c\C 2 /c q we could in 
principle have c\ — 0 or C 2 = 0. However, C 2 = 0 would mean that the gauge bosons would 
not couple at tree level to the Higgs, a scenario that is ruled out by experiment as shown 
in the previous section. Setting c\ — 0 we get, in Type II, ap = = 0 and 

4 = 44 bl = s 2 2 /t% bl = bl = t}sl C 2 = s}cl (4.1) 

In the left panel of figure 4 we show bjj — bL as a function of an — cll for Type II and 
a center of mass energy of 13 TeV with all rates at 10% (blue/black), 5% (red/dark-grey), 
and 1 % (cyan/light-grey) (in order to avoid the dependence on the phase conventions in 
choosing the range for the angles c^, we plot sgn(C) (sgn(C) bi) instead of ( bi ) with 
i = U,D,L). It is quite interesting to note that this scenario is still possible with the rates 
at 5% of the SM value at the LHC at 13 TeV. We have checked that this is still true at 
2 % and only when the accuracy reaches 1 % are we able to exclude the scenario. So far we 
have discussed ap = 0. Another interesting point is that when \cld\ —» 0, \b]j\ —> 1. The 
requirement that | ~ 1 implies that the couplings of the up-type quarks to the lightest 
Higgs take the form 

4 = a - 4) = (i -1/4)> 4 = 4 = 1/4 (4.2) 

while the coupling to massive gauge bosons is now 

c 2 = (ft -1 )/ft +1) = (1 - 4/(1 + 4 • 


(4-3) 




Type II 


Type II 


Q 

.a 



Figure 4. Left: sgn(C) bp = sgn(C) &£ as a function of sgn(C) au = sgn(C) cll for Type II and 
a center of mass energy of 13 TeV with all rates at 10% (blue/black), 5% (red/dark-grey), and 1% 
(cyan/light-grey). Right: same, but for sgn(C) by as a function of sgn(C) ajj . 


In the right panel of figure 4 we now show bjj as a function of ajj for Type II with the same 
colour code. We conclude from the plot that the constraint on the values of (a*/, bjj ) are 
already quite strong and will be much stronger in the future just taking into account the 
measurement of the rates. 

We would like to understand why ajj ~ 1 and by ~ 0, while the bounds are much 
looser for ay. We start by noting that the couplings impose different constraints on the up 
and down sectors. Indeed, from table 1 and eq. (2.9), 


C — si au 

Ru = -; Ri 2 = spau', R-i3 = tan f3 bjj , 

cp 

for the up sector, while 

C — Cp CL]J Cj3 

Rn = cr op; i?i 2 =-; R13 =- 

sp S(3 


for the down sector. In the first case, 


R \ 2 + R 2 


13 


2 

a u + -J < 

c p 


< 1 leads to 

1 


(4.4) 


(4.5) 


(4.6) 


Noting that the tan/3 > 1 constraint forces c# < l/y/2 and sp > 1/y/ 2, we find bjj < 1, 
while au < y/2. This is what we see in the right panel of figure 5, where in cyan we 
show points which are subject only to the theoretical constraints. We see that all points 
lye inside the ellipse in eq. (4.6). The constraint from the fiyy bound (orange/dark-grey 
points in the right panel of figure 5) then places the points on a section of that ellipse close 
to (ay, by) ~ (0,1). 

The situation is completely different for the down sector. Indeed, a similar analysis 
starting from eqs. (4.5) and R\ x + R 2 3 < 1, would lead to 


2 bjj 1 

a D + ~2 <U' 

° P C 0 


(4.7) 
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Figure 5. Left: Simulation points on the sgn (C) an versus sgn (C)bn plane. In cyan/light-grey 
(orange/dark-grey; blue/black) we show the points which pass all theoretical constraints (pass, in 
addition, the restriction from fiyy at 10%; pass, in addition, the restriction from all /i’s at 10%). 
Right: Same constraints on the sgn(C) au versus sgn(C) bu plane. 

Since cp can be very small, this entails no constraint at all, agreeing with the fact that 
the cyan/light-grey points in the left panel of figure 5 have no restriction. In contrast, it 
is the bound on fiyy which constrains the parameter space to the orange/dark-grey circle 
centered at (0,0). But now, the whole circle is allowed. 

The constraints from fiyy can be understood with simple arguments as follows. It 
was shown in [48], in the real 2HDM, that the limits on fiyy impose rather non trivial 
constraints on the coupling to fermions which, however, can be understood from simple 
trigonometry. Following the spirit of that article, we assume that the production is mainly 
due to gg -A h\ with an intermediate top in the triangle loop, and that the scalar decay 
width is dominated by the decay into bb. As a result, 

C 2 

Hvv ~ (a-u + l-Sfrjy) 2 _i_ ’ ( 4 -8) 

a D + °D 

where the approximate factor of 1.5 is what one would obtain either from a naive one- 
loop calculation 3 , or from a full HIGLU simulation [44]. Applying this formula, we obtain 
figure 6, where we have taken fiyy within 20% (orange/light-grey) or 10% (blue/black) 
of the SM, letting the angles vary freely within their theoretically allowed ranges. The 
similarity between the left (right) panes of figures 5 (a full model simulation) and figure 6 
(a simple trigonometric exercise) is uncanny. 

Further constraints are brought about by a second simple geometrical argument. They 
place all solutions close to (a, b) ~ (1,0) when C is close to unity. We use eqs. (4.4) to 
derive 

l = i?? 1 + J R ? 2 + i?? 3 = T^ 


3 See, for instances, equations (A7)-(A9) and (A14) in [6]. 




+ s 2 p a\j + tan z (3b(j, 


(4.9) 
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Figure 6. Left: sgn(C) bp = sgn(C) &£ as a function of sgn(C) ap = sgn(C) cll for Type II and 
a center of mass energy of 13 TeV with fiyy within 20% (orange/light-grey) and 10% (blue/black) 
of the SM. 

leading to 

+ >!(, = (4.10) 

This is a circle centered at ((7,0), which excludes most cyan/light-grey points on the right 
panel of figure 5. Since C is close to unity, and appears divided by tan/3 (which must be 
larger than one), the radius is almost zero, forcing ajj to lie close to C ~ 1, and bjj close to 
0. Including all channels at 10% restricts the region even further, as seen in the blue/black 
points on the right panel of figure 5. It is true that an equation similar to eq. (4.10) can be 
found for the down sector: 


(ap — C) 2 + b 2 D = tan 2 /3 [1 — C 2 ]. (4.11) 

However, the different placement of tan/3 is crucial. For intermediate to large tan/3, the 
tan 2 /3 factor in eq. (4.11) enhances the radius with respect to that allowed by the cot 2 /3 
factor in eq. (4.10). This explains the difference between the red/dark-grey points on the 
two panels in figure 4. 

We now turn to the constraints on the sina^-tan^ plane. When we choose fiyy > 0.9 
in the exact limit (|od|,|&d|) = (0,1), we obtain, using the approximation in eq. (3.5) 
tan/3 > 4.4. Because we are not in the exact limit, the bound we present in the left plot 
of figure 7 for tan/3 is closer to 3. The left panel of figure 7 shows tan/3 as a function of 
sina^ for Type II and a center of mass energy of 13 TeV with all rates at 10% (blue/black). 
In red/dark-grey we show the points with \ajj\ < 0.1 and ||&d| — 1| < 0.1 and in green 
IM < 0.05 and ||od| — 1| < 0.05. In the right panel, tan/3 is replaced by cosoq. These two 
plots allow us to distinguish the main features of the SM-like scenario, where (|ao|, |&d|) ~ 
(1,0) from the pseudoscalar scenario where (|od|, \bn\) ~ (0,1). In the SM-like scenario 
sina^ ~ 0, tan/3 is not constrained and the allowed values of sino^ grow with increasing 
cosoq. In the pseudoscalar scenario cosoq ~ 0, sino^ and tan/3 are strongly correlated 
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Figure 7. Left: tan/3 as a function of sirm 2 for Type II and a center of mass energy of 13 TeV 
with all rates at 10% (blue/black). In red/dark-grey we show the points with \au\ < 0.1 and 
\\b]j\ — 1| < 0.1 and in green \bu>\ < 0.05 and \\au>\ — 1| < 0.05 Right: same, with tan/3 replaced by 

COS OL\. 


and tan/3 has to be above ~ 3. Clearly, all values of ap and bp are allowed provided 
djj T bjj ~ 1. 



Figure 8. Left: sgn(C) bu as a function of sgn(C) ajj for Type I and a center of mass energy of 
13 TeV with all rates at 10% (blue/black) and 5% (red/dark-grey). Right: sgn(C) bu as a function 
of sgn(C) cll for LS and a center of mass energy of 13 TeV with all rates at 10% (blue/black) and 
5% (red/dark-grey). 

In the left panel of figure 8 we show bjj as a function of au for Type I and a center of 
mass energy of 13 TeV with all rates at 10% (blue/black) and 5% (red/dark-grey). In Type 
I this plot is valid for all Yukawa couplings, because au — cld — <Tl and bu = bu — bu- 
It is interesting that even at 10% there are points close to (a, b) — (0.5, 0.6) still allowed 
and no dramatic changes happen when we move to 5%. In the right plot we show 6 ^ as a 
function of au for LS with the same colour code. Here again the ( au , bu) = (0,1) scenario 
is still allowed both with 10% and 5% accuracy. However, as was previously shown, the 
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wrong sign limit is not allowed for the LS model [48, 49]. Nevertheless, in the C2HDM, the 
scalar component sgn (C) cll can reach values close to —0.8. Finally, for the up-type and 
down-type quarks, the plots are very similar to the one in the right panel of figure 4 for 
Type II. 

4.1 Direct measurements of the CP-violating angle 

Although precision measurements already constrain both the scalar and pseudoscalar com¬ 
ponents of the Yukawa couplings in the C2HDM, there is always the need for a direct (and 
thus, more model independent) measurement of the relative size of pseudoscalar to scalar 
components of the Yukawa couplings. The angle that measures this relative strength, fa, 
defined as 

tan fa — bijai i = U, D, L , (4-12) 

could in principle be measured for all Yukawa couplings. The experimental collaborations 
at CERN will certainly tackle this problem when the high luminosity stage is reached, 
through any variables able to measure the ratio of the pseudoscalar to scalar component of 
the Yukawa couplings. There are several proposals for a direct measurement of this ratio, 
which focus mainly on the tth and on the r + r~h couplings. Measurement of bjj/au were 
first proposed for pp tth in [50] and more recently reviewed in [51-53]. A proposal to 
probe the same vertex through the process pp hjj [54] was put forward in [55] and again 
more recently in [56]. In reference [56] an exclusion of fa > 40° ( fa > 25°) for a luminosity 
of 50 fb -1 (300 fb -1 ) was obtained for 14 TeV and assuming fa = 0 as the null hypothesis. 
A study of the r + r~h vertex was proposed in [57] and a detailed study taking into account 
the main backgrounds [58] lead to an estimate in the precision of Afa of 27° (14.3°) for a 
luminosity of 150 fb -1 (500 fb -1 ) and a center of mass energy of 14 Tev. 

Since in the C2HDM the couplings are not universal, one would need in principle three 
independents measurements, one for up-type quarks, one for down-type quarks and one 
for leptons. The number of independent measurements is of course model dependent. For 
Type I one such measurement is enough because the Yukawa couplings are universal. For 
all other Yukawa types we need two independent measurements. It is interesting to note 
that for model F, since the leptons and up-type quarks coupling to the Higgs are the same, 
a direct measurement of the hbb vertex is needed to probe the model. On the other hand, 
and again using model F as an example, a different result for fa and </> T would exclude 
model F (and also Type I). 

Let us first discuss what we can already say about the allowed range for the (j)jj = fa 
angle and what to expect by the end of the LHC’s run 2 using only the rates’ measurements. 
In figure 9 we show on the top row tan/3 (left), sin<a 2 (middle) and cosoq (right) as a 
function of (j)jj = tan ^{bjj/au) f° r Type I, with rates within 20% (green) and 5% (red/dark- 
grey) of the SM prediction. In the bottom row we present the same plots but for Type 
II. The green points are a good approximation for the allowed region after run 1, while 
the red/dark-grey points are a good prediction for the allowed space with the run 2 high 
luminosity results. The most striking features of the plots are the following. For Type I the 
angle <fiu = </>£> = </>£ is between —75° and 75° and this interval will be reduced to roughly 
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Figure 9. Top: tan/3 (left), sin 0:2 (middle) and cosaa (right) as a function of fa — tan 1 (bu/au ) 
for Type I with rates at 20% (green) and 5% (red/dark-grey). Bottom: same but for Type II. 


—45° and 45° provided the measured rates are in agreement with the SM predictions. For 
Type II only (j)jj is constrained; we get \(j)jj\ < 30° and the prediction of roughly \fa\ < 15° 
when rates are within 5% of the SM predictions. Since the Higgs couplings to top quarks 
are the same for all models, the angle that relates scalar and pseudoscalar components for 
this vertex is related to the lightest Higgs CP-violating angle by 

tan0t = — cp/si tan «2 => tan<a 2 = — si/cp tan fa . (4-13) 

The parameter space is restricted in such a way that high tan/3 implies low « 2 . Since s\ 
cannot be too small, it is clear from equation (4.13) that large tan/3 necessarily implies a 
small fa. This is clearly seen in the left top and bottom plots of figure 9 where for large tan /3 
the pseudoscalar component of the up-type quarks Yukawa coupling is very close to zero. 
Interestingly, for both Type I and Type II the values of tan /3 ~ 0(1) are the ones for which 
the angle fa is less constrained. These are exactly the values for which the the coupling tth 
has a maximum value (already considering the remaining constraints that disallow values 
of tan/3 below 1). Therefore, a direct measurement of fa could still be competitive with 
the rates measurement in Type I. 

Let us now move to the Yukawa versions that can have a zero scalar component not 
only at the end of run 1, but also at the end of run 2, if only the rates are considered. For 
definiteness we focus on Type II. As previously discussed, a direct measurement involving 
the vertex hr + r~ [57, 58] could lead to a precision in the measurement of fa , Afa , of 27° 
(14.3°) for a luminosity of 150 fb -1 (500 fb -1 ) and a center of mass energy of 14 Tev. In 
figure 10 (left) we show cosoq as a function of tan _1 (&£)/a£)) for Type II and a center of 
mass energy of 13 TeV with all rates at 10% (blue/black). In red/dark-grey we show the 
points with \afa < 0.1 and ||&d| —1| < 0.1 and in green \bfa < 0.05 and \\afa —1| < 0.05. In 
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Figure 10 . Left: cosaa as a function of tan _1 (6^/a^) for Type II and a center of mass energy of 
13 TeV with all rates at 10% (blue/black). In red/dark-grey we show the points with \cld\ < 0.1 and 
\\bo\ ~ 1| < 0.1 and in green \bp\ < 0.05 and \\cld\ — 1| < 0.05. Right: same, with cosaa replaced 
by sgn(C) a D . 


the right panel cosoa is replaced by sgn(C) op. It is clear that the SM-like scenario sgn(C) 
(a/), bjj) = (1,0) is easily distinguishable from the (0,1) scenario. In fact, a measurement 
of </> T even if not very precise would easily exclude one of the scenarios. Obviously, all other 
scenarios in between these two will need more precision (and other measurements) to find 
the values of scalar and pseudoscalar components. The r + r~h angle is related to OL2 as 

tan</> T = —sp/ci tan c ^2 => tan<a 2 = —c\/sp tan </> T (4-14) 

and therefore a measurement of the angle </> T does not directly constrain the angle <a 2 . In 
fact, the measurement gives a relation between the three angles. A measurement of <pt and 
<j) T would give us two independent relations to determine the three angles. 

4.2 Constraints from EDM 

Models with CP violation are constrained by bounds on the electric dipole moments (EDMs) 
of neutrons, atoms and molecules. Recently the ACME Collaboration [59] improved the 
bounds on the electron EDM by looking at the EDM of the ThO molecule. This prompted 
several groups to look again at the subject. For what concerns us here, the complex 2HDM, 
several analyses have been performed recently [7, 10, 60-63]. In ref. [7] it was found that 
the most stringent limits are obtained from the ThO experiment, except in cases where 
there are cancellations among the neutral scalars. These cancellations were pointed out in 
[62, 63] and arise due to orthogonality of the R matrix in the case of almost degenerate 
scalars [9]. So far, there is no complete scan of EDM in the C2HDM; only some benchmark 
points have been considered, making it difficult to see when these cancellations are present. 

What can be learned from these studies is that the EDMs are very important and 
their effect in the C2HDM has to be taken in account in a systematic way, in the sense 
that, for each point in the scan, the EDMs have to be calculated and compared with the 
experimental bounds. However, for the purpose of the studies in this work and for the 
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present experimental sensitivity, this is not required. This is because we are looking at 
scenarios where the couplings of the up-type sector (top quark) are very close to the SM 
and the differences, still allowed by the LHC data, are in the couplings of the down-type 
sector; the tau lepton and bottom quark. As was shown in ref. [10], while the pseudoscalar 
coupling of the top quark is very much constrained (in our notation \bjj\ < 0.01), the 
corresponding couplings for the b quark and tau lepton are less constrained by the EDMs 
than by the LHC data. Since we are taking in account the collider data, our scenarios are 
in agreement with the present experimental data. But one should keep in mind that, as 
pointed out in refs. [10, 64], the future bounds from the EDMs can alter this situation. In 
the future, the interplay between the EDM bounds and the data from the LHC Run 2 will 
pose relevant new constraints in the complex 2HDM in general, and in particular for the 
scenarios presented in this work. 

5 Conclusions 

We discuss the present status of the allowed parameter space of the complex two-Higgs 
doublet model where we have considered all pre-LHC plus the theoretical constraints on 
the model. We have also taken into account the bounds arising from assuming that the 
lightest scalar of the model is 125 GeV Higgs boson discovered at the LHC. We have shown 
that the parameter space is already quite constrained and recovered all the limits on the 
couplings of a CP-conserving 125 GeV Higgs. The allowed space for some variables, as for 
example for the tan /? parameter, is now increased as a natural consequence of having a 
larger number of variables to fit the data as compared to the CP-conserving case. 

The core of the work is the discussion of scenarios where the scalar component of the 
Yukawa couplings of the lightest Higgs to down-type quarks and/or to leptons can vanish. 
In these scenarios, that can occur for Type II, F and LS, the pseudoscalar component 
plays the role of the scalar component in assuring the measured rates at the LHC. A direct 
measurement of the angle that gauges the ratio of pseudoscalar to scalar components in 
the tth vertex, </>£, will probably help to further constrain this ratio. However, it is the 
measurement of </> T , the angle for the r + r~h vertex, that will allow to rule out the scenario 
of a vanishing scalar even with a poor accuracy. We have also noted that for the F model 
only a direct measurement of </>£> in a process involving the bbh vertex would be able to 
probe the vanishing scalar scenario. Finally a future linear collider [65, 66] will certainly 
help to further probe the vanishing scalar scenarios. 
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